Nonlinear Dynamics of a Bose-Einstein Condensate in a Magnetic Waveguide 
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We have studied the internal and external dynamics of a Bose-Einstein condensate in an anhar- 
monic magnetic waveguide. An oscillating condensate experiences a strong coupling between the 
center of mass motion and the internal collective modes. Due to the anharmonicity of the magnetic 
potential, not only the center of mass motion shows harmonic frequency generation, but also the 
internal dynamics exhibit nonlinear frequency mixing. Thereby, the condensate shows shape oscil- 
lations with an extremely large change in the aspect ratio of up to a factor of 10. We describe the 
data with a theoretical model to high accuracy. For strong excitations we test the experimental 
data for indications of a chaotic behavior. 

PACS numbers: 03.75.Kk, 03.75.Be, 05.45.-a 



Since the first realization of Bose-Einstein condensates 
in 1995 ideas have been discussed how to exploit this new 
form of coherent matter for atom optics. Atom inter fer- 
ometry with single atoms has been extensively investi- 
gated in the past decade [l[ and finds important applica- 
tions for high precision measurements of gravitational or 
rotational forces (2, La] • It was also possible to realize and 
demonstrate a number of atom optical elements such as 
mirrors Q , gratings , beamsplitters Q and waveguides 
7,8]. It is a fascinating vision to integrate these elements 
on a microfabricated atom chip. Extremely sensitive 
chip based sensors for forces and accelerations as well as 
for applications in quantum computing are conceivable. 
Since Bose-Einstein condensates have been loaded into 
microfabricated traps 0, this vision has changed into 
a very real research topic. Integrated atom optics with 
condensates combine the nonlinear interaction of a Bose- 
Einstein condensates with trapping potentials, whose ge- 
ometry is in general anharmonic. This is in contrast 
to both single atom experiments, where the interaction 
plays no important role, and conventional experimental 
arrangements for condensates, where the trapping poten- 
tial is typically harmonic and consequently the center of 
mass (CM) motion is decoupled from the internal dynam- 
ics |ll|. In microtraps one therefore expects not only a 
coupling between the external and internal dynamics but 
also pronounced nonlinear effects including even chaotic 
dynamics for large excitation amplitudes. The study of 
the dynamics of a condensate in micropotentials is an 
extension of previous work on collective excitations and 
provides essential information for future applications for 
integrated atom optics. 

We present a detailed experimental study of the con- 
densate dynamics in the most fundamental atom optical 
element — a magnetic waveguide. We regard waveguides 
which are closed at both ends such that the generic mo- 
tion of the condensate is a periodic oscillation. The linear 
propagation of Bose-Einstein condensates along a non 
terminated guide was recently investigated in Refs. [l^ 



and We have measured the time evolution of the 

CM motion and of the aspect ratio i. e. the ratio between 
the radial and the axial condensate radii. In contrast to 
previous experimental work on collective modes of a con- 
densate, we realize a closed system with conserved en- 
ergy and a self-sustained excitation scheme. This system 
allows us to realize high amplitude oscillations of the as- 
pect ratio which exceed previously reported values by far. 
To describe the experimental data we have developed a 
theoretical model which is described in Ref. |15j . The 
model is based on the separation of the CM motion and 
reduces the internal dynamics to that of a condensate 
in a fictitious time-dependent harmonic potential. Such 
a scenario is known to be exactly solvable in Thomas- 
Fermi approximation 0, . The model also allows us 
to identify an experimentally accessible regime where the 
dynamics become chaotic. 

The waveguide potential is generated with microfab- 
ricated copper conductors at the surface of a chip. The 
chip is mounted upside down in a vacuum chamber at a 
background pressure of 10 -11 mbar. The 87 Rb atoms are 
collected in a six-beam magneto-optical trap, pumped 
into the (F = 2, mp = 2) ground state and trapped in a 
spherical magnetic quadrupole trap. The atoms are then 
adiabatically transferred into the waveguide potential at 
the chip surface by gradually changing the magnetic field 
geometry 0, [HI • By forced evaporation, Bose-Einstein 
condensates with 100.000 atoms are produced. For the 
presented experiments, the trapping frequencies are ad- 
justed to 110 Hz in the radial and 8 Hz in the axial di- 
rection. The shape of the axial potential is sketched in 
Fig.^i. The waveguide potential is characterized by an 
almost constant radial oscillation frequency (the maximal 
deviations are smaller than 10 percent) and a strongly 
varying axial curvature. With a numerical calculation of 
the magnetic field, we determine the axial potential U (x) 
to 
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FIG. 2: Evolution of the aspect ratio of the condensate for 
the data series A (a), B (b), and C (c) (dots) and theoretical 
model (solid line), (d) Plot (c) with a logarithmic scale, (e) 
All data points of series B. 



FIG. 1: Oscillation of the condensate in an anharmonic 
waveguide. The condensate is trapped below the microstruc- 
ture which generates the magnetic field for the trapping po- 
tential, (a) Magnetic field in the trap minimum along the ax- 
ial direction. The black dots indicate the starting positions of 
the oscillation for each experimental series. The circle marks 
the initial position without displacement, (b) Absorption im- 
ages of the condensate after 20 ms time of flight (series C, 
10ms intervals). Due to the time of flight, the amplitude of 
the oscillation is enhanced with respect to oscillation in the 
trap (also a phase shift occurs). For condensates marked with 
an asterisk the absorption images show diffraction and the as- 
pect ratio can not be determined properly, (c) CM motion 
of the condensate (series B): experimental data (dots) and fit 
with three sinusoidal functions (line). 



atom number. To exclude influences of the nearby cop- 
per conductors [HI E(j , the experiment was carried out 
in a distance to the surface of 220 fim. 

Our theoretical model describes the condensate with 
a set of two coupled differential equation systems. The 
first system are the classical equations of motion for the 
CM position of the condensate 



raR= -W(R), 



(2) 



where m is the mass of the atom, R is the CM position 
and U is the confining potential. In Thomas-Fermi ap- 
proximation, the density distribution of a condensate in 
a time-dependent harmonic potential can be written as 
follows 



where x has to be quoted in millimeter and /i is the mag- 
netic moment of the atom. To induce the oscillation 
we displace the condensate within 400 ms with a mag- 
netic field gradient and switch off the field gradient non- 
adiabatically (< 60 /is). The experiment was performed 
for three different initial displacements d: 0.45 mm (se- 
ries A, duration t& = 1.015s), 0.63mm (B, t& — 2.7s) 
and 0.79mm (C, t c = 0.765s). In FigJ^b the first 1.5 
oscillation periods of series C are shown. To determine 
the CM position and the aspect ratio, the absorption im- 
ages are fitted to a Thomas- Fermi profile. The aspect 
ratio bears the advantage of being independent from the 



n(x, t) 



g X 1 (t)X 2 (t)X 3 (t) 



(3) 

where fi is the chemical potential and g is the coupling 
constant. The Xi(t) contain the complete time evolution 
of the Thomas-Fermi radii ri(t) = r^Xiit). The are 
governed by three ordinary differential equations 
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Here, U0i(t) are the time-dependent oscillation frequen- 
cies of the harmonic potential. If the modulation of the 
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potential starts at t = 0, the initial conditions for the pa- 
rameters are a^(0) = ujiq and 7*^(0) = with A^(0) = 1. 
The anharmonicity of the potential now provides a cou- 
pling of the CM motion and the time evolution of the 
Thomas-Fermi radii, which is given by 

^ {t) = i^ u{K) ' i = 1 ' 2 ' 3 - ( 5 ) 

To compare the experimental data with the model, (0), 
J3J) and (0 are numerically solved with the initial condi- 
tions Ai(0) = 1 and R(0) = (d,0,0). 

The experimental data for the CM motion of series B 
are presented in Fig.^ . We cannot extract any damp- 
ing and the data are consistent with a (1/e) -lifetime of 
at least 20 minutes. Due to the anharmonicity of the 
potential (QJ, the frequency spectrum of the CM mo- 
tion shows second and third harmonic generation. For 
series C with the largest initial displacement, the non- 
linear effects are most pronounced and lead to a relative 
amplitude of 5 percent for the second harmonic and 1 
percent for the third harmonic frequency with respect 
to the amplitude for the fundamental frequency For 
small amplitudes amounts to 8.5 Hz but is reduced 
with increasing displacement (A, B, and C) by 4, 6, and 
11 percent respectively. Because of the extremely low 
damping, the CM motion is suited to detect small forces, 
which may shift the trapping potential. An estimate for 
our experimental parameters indicates a resulting sensi- 
tivity of F/m — 10 -4 g, with the gravitational accelera- 
tion g. 

In Fig.[2M21i the measured aspect ratio is compared 
with the theoretical model. While the theory contains 
no free parameter, a constant multiplication factor of 1.2 
is needed to describe the data accurately. This factor 
arises from the fact, that the condensate is radially com- 
pressed during the switch-off of the magnetic trap which 
- due to the mean field repulsion — leads to a larger 
radial size in the time-of- flight image [2]]. The agree- 
ment between experiment and theory is remarkable and 
proves true the picture that the external motion in an an- 
harmonic potential leads to a time-dependent harmonic 
potential in the rest frame. In Fig.^D, the condensates 
marked with an asterisk cannot be spatially resolved by 
our absorption imaging system, indicating an axial ex- 
tension of the condensate of less than 15 jam. Therefore, 
the condensate may in fact follow the strong dynamics, 
which is predicted by the theory (FigEJi). For the present 
data, we measure a change in the aspect ratio by more 
than a factor of 10, the largest value reported so far. The 
data of series B which cover 2.7 s are promising to find 
effects arising from loss of atoms (FigEt). At the end of 
the series, the atom number has dropped by more than 
a factor of 15 with only 3000 atoms left in the conden- 
sate. The systematic increase of the aspect ratio is due 
to the crossover to the ID regime as defined in Ref. |22j . 
The Thomas-Fermi approximation breaks down because 
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FIG. 3: Frequency spectra of the aspect ratio. The frequen- 
cies of the labeled peaks are: (i) fundamental frequency z^n, 
(ii) second harmonic frequency 2z/n, (iii) quadrupole resonance 
frequency v q ~ yj 5 / 2 z^n , (iv) difference frequency 2v§ — u q , (v) 
sum frequency 2vo + v q and (vi) fourth harmonic frequency 
4z/ . 

the chemical potential of the condensate becomes smaller 
than the radial kinetic energy, which then dominates the 
radial expansion of the condensate. 

We have performed a frequency analysis of the aspect 
ratio with an algorithm suitable for unequally spaced or 
missing data 23]. The results are shown in Fig.03 The 
spectra show three major frequencies: the fundamental 
frequency (peak (i) in Fig.|3|), the second harmonic 
(ii) and the resonance frequency v q of the lowest lying 
collective mode (iii). The latter does not deviate more 
than 4 percent from its small amplitude value for strongly 
anisotropic traps which is given by y/5 / 2 zv . Further- 
more, the spectrum for series C (Fig.E) shows several 
mixed frequencies. Two of them, the sum and the dif- 
ference between the second harmonic frequency and the 
lowest resonance frequency 2zv — v q and 2zv + v q have 
been indentified ((iv) and (v)). This gives evidence for 
the nonlinear coupling of the collective excitations. Other 
than the coupled modes observed by Hechenblaikner et 
at. |24|, the mode coupling in our experiment is non- 
resonant, and the mixed frequencies vanish in the limit 
of small amplitudes. Because the condensate performs a 
free oscillation in an external potential the total energy 
of the system is conserved. Summing up only the con- 
tributions for the internal energy of the condensate, one 
finds that the sum is not constant. Adding the kinetic 
and potential energy of the CM motion, the energy is 
conserved, which can be verified by using a Hamiltonian 
formalism Therefore, not only the internal dynam- 
ics is affected by the external motion, but there is also 
an energy exchange between the internal dynamics and 
the CM motion. In our experiment this effect leads to a 
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FIG. 4: (a) Oscillation for an initial displacement of 0.94 mm. 
The absorption images are taken after 20 ms time of flight. 
The vertical dotted line indicates the surface of a wire which 
narrows the display window, (b) Lyapunov exponent: An = 
as a function of the initial displacement 
d. do is the initial Euclidian distance between two neighboring 
trajectories, r is the duration of one iteration and N is the 
number of iterations. The Lyapunov exponent is defined in 
the limit A = limjv-»-oo An. Each pair of trajectories was 
numerically solved for 200 s (do = 0.1 /xm, r = 1ms, N = 
2 x 10 5 ). Above a displacement of 0.67 mm the numerics are 
no longer stable for the total propagation time. 



punov exponent becomes positive and grows linearly. As 
the total energy of the system is conserved, one can eas- 
ily test the numerical stability and exclude any numeri- 
cal artefacts. The corresponding Poincare maps and the 
frequency spectra of the trajectories indicate the same 
results. This gives strong evidence, that series C was 
performed in a regime where the longtime behavior of 
the trajectory is chaotic. 

In conclusion, we have studied the dynamics of a Bose- 
Einstein condensate oscillating in an anharmonic trap- 
ping potential. We have identified several nonlinear ef- 
fects including second and third harmonic generation of 
the center of mass motion, a coupling of the internal and 
external dynamics and nonlinear mode mixing. The the- 
oretical model describes these effects with high accuracy. 
The model predicts a chaotic dynamic above a critical 
oscillation amplitude and the experimental data give in- 
dication for a chaotic motion too. Considering applica- 
tions of magnetic microtraps and integrated atom op- 
tics, we have shown that the motion of a condensate in 
a magnetic waveguide necessarily leads to excitations of 
the inner degrees of freedom. A thorough understanding 
of these effects plays a key role to exploit the nonlinearity 
to build high precision sensors. It also shows fundamen- 
tal limits for the use of Bose-Einstein condensates for 
technical applications. 

We gratefully acknowledge financial support from 
the Deutsche Forschungsgemeinschaft under Grant No. 
Zi/419-5. 



theoretically calculated amplitude variation for the CM 
motion of 1/4000, which is too small to be detected. 

Fig.^i shows absorption images for an initial displace- 
ment of 0.94 mm. The data show a strong axial excita- 
tion of the condensate, which raises the question of the 
onset of chaos. A signature of a chaotic behavior in ex- 
perimental data is an enhanced noise in the frequency 
spectrum of the data. The spectrum of the aspect ratio 
for series C in Fig.[S| clearly shows an enhanced back- 
ground, compared to the two other series. However, this 
should not be regarded as an unchallengeable experimen- 
tal proof. Such a proof would require a longer data 
set and an enhanced resolution of the imaging system. 
Nevertheless, the theoretical model can be analyzed for 
the experimental parameters. One of the strongest evi- 
dences for a chaotic dynamic is a positive Lyapunov ex- 
ponent, which indicates an exponential sensitivity to the 
initial conditions. For this purpose neighboring trajecto- 
ries for slightly different initial displacements have been 
simulated over 200 s. Fig.Q] shows the evolution of the 
Lyapunov exponent with increasing displacement. For a 
displacement of more than 0.63 mm, the calculated Lya- 
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